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Abstract. The blue supergiant Sher25 in the massive Galactic cluster NGC3603 is surrounded by a striking 
emission line nebula. The nebula contains an equatorial ring and probable bi-polar outflows, and is similar in 
I morphology, mass and kinematics to the shell now visible around SN1987A. It has been suggested that both 

1^-^ . nebulae were ejected while Sher25 and the progenitor of SN1987A were in previous red supergiant phases. In the 

' case of Sher 25 this is based on the qualitative strengths of nebular [N ll] emission which is indicative of nitrogen 

fS) , enriched gas. This gas may have been dredged up to the stellar surface by convective mixing during a previous 

■ red supergiant phase. We present optical high- resolution spectra of Sher 25 and a model photosphere and unified 

stellar wind analysis which determines the atmospheric parameters, mass- loss rate and photospheric abundances 
for C, N, O, Mg, and Si. We compare these results, in particular CNO, to other Galactic B-type supergiants and 
find that Sher 25 does not appear extreme or abnormal in terms of its photospheric nitrogen abundance. The C/N 
O ' and N/O ratios are compared to surface abundances predicted by stellar evolutionary calculations which assume 

the star has gone through a red supergiant phase and convective dredge-up. In particular we find that the N/O 
^ ' abundance is incompatible with the star having a previous red-supergiant phase, and that the nebulae is likely to 

have been ejected while the star was a blue supergiant. The results are compatible with some degree of rotationally 
induced mixing having occurred while the star was on or near the main- sequence. This is similar to what has 
recently been found for nebulae surrounding LBVs. In addition our wind analysis suggests the star currently has a 
, . , relatively normal mass-loss rate in comparison with other Galactic B-type supergiants and sits comfortably within 

' the wind momentum-luminosity relationship. In light of the evidence regarding massive evolved early-type stars 

in the Galaxy we suggest there is no object which shows any evidence of having had a previous red supergiant 
phase and hence of undergoing blue loops in the HR diagram. 

Key words, stars: supernovae — stars: individual(Sher25) — stars: supergiants — stars: abundances — stars: 
winds, outflows 



1. Introduction 

The study of hot luminous B-type supergiant stars is of 
fundamental importance for understanding the late evo- 
lutionary stages of massive stars and the chemical evolu- 
tion of our own and external galaxies. Massive stars enrich 
the interstellar medium with the heavy chemical elements 
during all stages of their evolution, from fast winds in 
their early hot stages through cool dense outflows from 
red supergiants and finally and most spectacularly during 
their deaths in core-collapse supernovae explosions. The 
fact that the progenitor of SN1987A was identified as a 
hot B-type supergiant makes them extremely interesting 
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objects, being the final stage in the life of at least some 
very massive stars. 

NGC 3603 is a massive H ii region ionized by a young, 
luminous cluster located approximately 7kpc from the Sun 
in the fourth galactic quadrant (/=291°.62, 6=-0°.52). It 
has been described as the most massive optically visible 
giant Hil region in the Galaxy (Goss & Radhakrishnan 
1969). Given its relative proximity and the large number 
of massive stars it contains, this cluster provides an op- 
portunity to study the upper mass region of the IMF and 
the evolutionary processes of very massive stars. It con- 
tains approximately 50 O and early B-type stars and the 
age of the original starburst is estimated somewhere be- 
tween 2-4 Myr (Melnick et al. |1989|, Moffat |l 9831 Crowther 



2 



Smartt et al.: The blue supergiant Sher25 



& Dessart 1998). The cluster includes the trapezium- type 
system HD97950 which has three WR stars identified and 
six 03 stars all within 0.3 pc (Drissen et al. 1995). Of 
particular interest is the star no. 25 reported first by 
Sher (1965). Moffat (1982) presented a moderate resolu- 
tion spectrum of this blue supergiant and classified it as 
a B1.5Iab and Van den Bergh (1978) has described it as 
the most highly evolved star in the cluster. 

The main reason for the interest in this star is the 
spectacular circumstellar nebula, discovered by Brandner 



et al. (1997a) and beautifully illustrated in HST press 
release images STScI-PRC99-20^ It consists of a ring- 
shaped equatorial emission line nebula (expansion velocity 
of 20 km s^^) and probable bipolar outflows to the north- 
east and to the southwest (expansion velocity of 83 km 
The Ha+ [N ii] and NIR images of Brandner et al. 



( |19974 |1997b| ) indicate that both the ring and the bi- 
polar ejecta have roughly the same dynamical age which 
would suggest a common origin in a stellar outburst ap- 
proximately 6600 yrs ago. The [Nil] /Ha emission line ra- 
tio of the ring appeared to be higher than that of the 
background Hii region, which Brandner et al. have qual- 
itatively interpreted as evidence for nitrogen enrichment 
in the ejecta and hence possibly that the gas has been 
processed through the CN and ON-cycles in the core of 
Sher 25 during its main-sequence phase. However there 
are no quantitative determinations of electron tempera- 
ture, density and hence abundances in either the ring or 
bi-polar lobes. Sher 25 is currently a hot blue supergiant 
and Brandner et al. have suggested that the star lived 
briefly for a period as a red supergiant during which time 
it ejected the nebula. In this picture Sher 25 has under- 
gone a "blue-loop" in the HR diagram, flrst evolving to 
the red after main-sequence core-H exhaustion and then 
back again to the B-type supergiant region. However these 
conclusions are based mainly on the existence and mor- 
phology of the nebula and the suggestion that it may be 
enriched with nitrogen. How massive stars generally evolve 
in the upper regions of the HR diagram is still not well un- 
derstood. Theoretical models have many free parameters 
which can affect the path of a star after (and in some cases 
during) the main-sequence lifetime, such as mass-loss, ro- 
tation, convection and overshooting treatment, and metal- 
licity. Observational evidence on the evolutionary status 
and history of evolved massive stars such as blue and red 
supergiants, LBVs and WR stars provide necessary con- 
straints on the models. Hence determining the nature of 
the circumstellar gas and the evolutionary phase of Sher 25 
is of great interest in this context. In addition it has been 
suggested (by Brandner et al. 1997b| ) that Sher 25 could 
be a twin of the progenitor of SN 1987A, which was a BSIa 
supergiant (Walborn et al. 1989 ). This is due to the exis- 
tence of a similar ring structure now seen around SN1987A 
which in existence before the SN explosion and hence was 
ejected by the progenitor star. Fransson et al. ( |1989D sug- 
gested that it had passed through a previous RSG phase 



given the CNO line strengths in the early-time UV spec- 
tra, although the uncertainties are large. Therefore the 
star Sher 25 could be a possible supernova Type II pro- 
genitor. 

In the wider picture the existence and relative num- 
ber of blue supergiants in different mass ranges needs to 
be understood to provide stronger constraints on evolu- 
tionary theory. The problem with predicting the correct 
number ratio of blue-red supergiants in different metal- 
licity regimes was first pointed out by Langer & Maeder 



(1995). The B/R-ratio appears to increase with metallic- 
ity and such a trend has not been gene rally reproduced 
by theory. Recently Maeder & Meynet ( ^001 ) have sug- 
gested that introducing rotation into the calculations does 
go some way to improving the comparison with observa- 
tions. If massive stars do go through blue- loops then the 
blue-supergiant population may be made up of pre- and 
po st red -supergiant objects (e.g. see models of Schaller et 
al. 1992). However distinguishing between two populations 
and determining their ratios has been a notoriously diffi- 
cult problem. One method that has been tried is using stel- 
lar surface chemistry as an indicator of evolutionary sta- 
tus. Qualitatively for example Jaschek & Jaschek (1967) 
and later Walborn ( 1972| ) suggested that the carbon and 
nitrogen line strengths in OB stars could be used to define 
separate populations, which led to the OBN/OBC spec- 
tral classification. In this scheme the OBN stars generally 
have stronger nitrogen and weaker carbon lines, and vice- 
versa for the OBC stars. Also there is a wide variation in 
the strength of the differences found between stars. These 
stars are usually explained by invoking some mechanism 
that will mix core processed gas to the surface. This could 
be either dredge up by convective mixing during a previous 
red supergiant phase, or alternatively rotationally induced 
mixing during and just after the main-sequence phase 
could provide an efficient mechanism (e.g. Denissenknov 



1994, Talon et al. 1997, Maeder & Meynet 2001) 



^ http://oposite.stsci.edu/pubinfo/PR/1999/20/index.html 



As the nebulae around Sher 25 could be suggestive of it 
having definitely gone through a previous red-supergiant 
phase, it is of great interest to study its atmospheric and 
wind parameters in detail and deduce its surface chemi- 
cal composition. These can then be compared with values 
estimated from different evolutionary scenarios, and also 
compared to the general population of Galactic B-type 
supergiants. A great deal of effort has recently been in- 
vested in understanding Galactic B-type supergiants and 
two papers in particular have set cornerstones in the field. 
The first is the non-LTE model atmosphere study of atmo- 
spheric parameters, line strengths and abundance trends 
of McErlean et al. ( 1999| , hereafter Paper I), which stud- 
ied a large sample of B-type supergiants in a homogeneous 
way. The difficulties in analysis methods and ways of mit- 
igating errors were presented which provides an excellent 
starting point for future detailed analysis of single objects 
or homogeneous groups. The second is the extensive study 
of the wind parameters o f four teen early B-type super- 
giants by Kudritzki et al. ( 1999 hereafter Paper II). This 
used non-LTE unified model atmospheres to investigate 
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the Ha line profiles and determine stellar wind properties 
for a subset of the stars in Paper I (and some A-type su- 
pergiants in addition). A tight relationship between the 
stellar wind momentum and luminosity was found, and 
the Balmcr line profiles were quantitatively reproduced. 

In this paper we present a non-LTE model photosphere 
analysis to deduce the surface parameters and abundances 
of Sher25. In addition, a non-LTE unified model atmo- 
sphere is used to analyse the Ha emission line and char- 
acterise the stellar wind. We place the results of these 
analysis in the general context of the previously analysed 
large sets of Galactic B-type supergiants. 



2. Observational data 

Spectra of Sher25 {V = 12.3,5 ^ 13.7; Brandner et al. 



1997a ) were taken with the UCLES echelle spectrometer 
at the Anglo-Australian Telescope on the night of 28th 
June 1999. The 79 Imm^i grating and the MIT2 CCD 
detector were used. The CCD was of format 4kx2k with 
15/im pixels, which results in a 1" entrance slit project- 
ing onto four pixels on the detector. To avoid this over- 
sampling and increase the signal-to-noise per resolution 
element, we binned the CCD on readout to 2kxlk. Two 
wavelength regions were observed, one centred on 4335A 
and the other on 6560A. The slit width was kept constant 
at 1", and the FWHM of the arc lines at the centres of 
these settings were measured at O.llA and 0.17A respec- 
tively. One exposure of 1200s was taken in the red setting, 
and 4x 1200s were taken in the blue; exposures of compar- 
ison arc lamps were taken immediately after each science 
set of exposures. The purpose of the red setting was mainly 
to get the Ha profile, while the blue covers the higher 
lines in the Balmer series, and all of the He i and weaker 
metal lines required for atmospheric diagnostics. The 2D 
frames were bias corrected, fiat-fielded, and the spectral 
orders were optimally extracted using standard routines 
within IRAF. After wavelength calibration the four blue 
frames were cross-correlated (taking the best S/N spec- 
trum as a reference) to check for velocity shifts. None were 
found at a significance greater than the estimated error of 



■^6 km 



so the spectra were co-added. The final blue 



spectra had a dispersion of 0.06Apix~^ in the central or- 
der, which considerably over-samples the metal line pro- 
file widths. The mean FWHM of several of the stronger 
metal lines (in the region 4550-4620A) was measured at 
I.SA (98 km s~^), hence the spectra were rebinned to a 
dispersion of 0.3Apix~^. At the furthest useful blue wave- 
length (3915A) this corresponds to 4 resolution elements 
sampling the FWHM which is adequate for line profile fit- 
ting and significantly increased the S/N in the continuum. 
The MIT2 CCD is too large for the window of the UCLES 
camera, hence the ends of each echelle order are severely 
vignetted. The ends of each order were thus trimmed, and 
the final useful spectrum has almost complete coverage of 
3915-5280A across 15 orders, with gaps of 5-I5A appear- 
ing only between the last 5. The signal- noise-ratio ranges 



between 50 at the edge orders to 70 in the central few 
orders. 

The spectra were normalised by fitting cubic splines to 
continuum regions free from absorption lines in iraf and 
afterwards were transferred to STA RLINK spectral analy- 
sis program DIPSO (Howarth et al. 1998 ). The equivalent 
widths of the metal lines and non-diffuse lines of neutral 
helium were measured by the non-linear least squares fit- 
ting of single or multiple Gaussian profiles to the nor- 
malised spectra. The hydrogen and diffuse helium lines 
were not measured in this manner, but the normalised pro- 
files were extracted directly for comparison with Galactic 
standards and with theoretical profiles. 

In order to compare the atmospheric parameters, pho- 
tospheric abundances and wind parameters with other 
typical B-type supergiants, three bright stars of similar 
spectral type (lying within the solar neighbourhood) were 
chosen as comparison standards. These were HD2905, 
HD 14956 and HD 13854, which have very similar He, 5,7 
lines and Siiv/Siiii line strength ratios as Sher25. These 
stars have been observed with the single-order spectro- 
graph ISIS on the WHT, and the spectra have a coverage 
of 3900-4735A covering the main diagnostic lines of in- 
terest (see Smartt et al. 2001 where the details of the in- 
strument setup is described). The S/N of these spectra are 
very high (in excess of 200) and the line profiles are well re- 
solved, hence they make excellent comparison templates. 
Equivalent widths for these stars were also measured in 
a similar method described above. These three stars were 



also part of the Lennon et al. (1993) atlas, which covers the 
region 3950-5000A, and the Ha region 6260-6870A. Line 
strength measurements in wavelength regions not covered 
by the WHT ISIS spectra (mainly 4735-500GA) were taken 
from this source. The final values of equivalent widths for 
each star can be found in Table 2. 

Table 2. The atmospheric parameters derived for Sher 25 
and the three comparison supergiants, with the homoge- 
neous method described in Sect.|[ 



Star 




logs 




V sin i 




(K) 


(cgs) 


(kms^^) 


(kms~^) 


Sher 25 


22 300 ±1000 


2.60 ±0.1 


15 ±5 


65 


HD2905 


23 500 ±1000 


2.70 ±0.1 


11 ±5 


80 


HD13854 


22 000 ±1000 


2.65 ±0.1 


13 ±5 


80 


HD 14956 


21250 ±1000 


2.55 ±0.1 


12 ±5 


75 



3. Model atmosphere and wind analysis 

3.1. Photospheric analysis and abundances 

The photosphere of Sher 25, and subsequently the three 
comparison stars, were modeled using the techniques de- 
scribed in det ail in Paper I and briefly summarized in 
Smartt et al. ( ^OOla ). A grid of non-LTE model atmo- 
spheres has previously been generated using the code 
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Table 1. Equivalent widths for metal lines of Sher25 and comparison stars given in mA units. 



Ion and Line 


HD2905 


HD13854 


HD 14956 


Sher 25 


C II 4267.02 


111 


155 


128 


125 


& 4267.27 










C II 6578.10 




135 


105 


110 


N II 3995.00 


101 


308 


443 


484 


N II 4035.08 


25 


89 


101 


95 


N II 4236.86 


28 


78 


63 


75 


& 4236.98 










N II 4241.18 


53 


119 


88 


165 


& 4241.78 










N II 4447.03 


70 


106 


172 


134 


N II 4607.16 


21 


134 


247 


265 


N II 4613.87 


16 


91 


188 


160 


N II 4630.54 


79 


301 


456 


398 


N II 4803.29 




122 


187 


107 


O II 4069.62 


249 


235 


171 


167 


& 4069.89 










O II 4072.16 


213 


220 


190 


122 


O II 4075.86 


235 


241 


204 


181 


O II 4132.80 


85 


81 


63 


65 


II 4317.14 


204 


192 


135 


177 


II 4319.63 


205 


210 


155 


210 


& 4319.93 










II 4345.56 


264 


202 


147 


184 


II 4347.42 




112 


77 


91 


O II 4349.43 


397 


347 


267 


297 


II 4351.26 


183 


127 


85 


103 


& 4351.50 










II 4366.89 


236 


217 


170 


235 


II 4414.90 


303 


338 


257 


291 


O II 4416.97 


240 


270 


211 


185 


O II 4590.97 


207 


173 


147 


170 


II 4595.96 


177 


137 


159 


134 


& 4596.18 










II 4638.86 


242 


212 


126 


159 


II 4661.63 


259 


233 


202 


194 


II 4673.74 


88 


57 


53 


60 


II 4676.24 


234 


222 


150 


253 


Mg II 4481.13 


170 


244 


233 


205 


& 4481.33 










Si IV 4116.10 


202 


141 


111 


218 


Si III 4552.62 


397 


455 


485 


555 


Si III 4568.28 


348 


390 


438 


437 


Si III 4575.22 


219 


241 


293 


266 


Si III 4813.30 


50 


65 


55 


49 


Si III 4819.72 


80 


90 


100 


96 


Si III 4829.96 




90 


85 


91 


S III 4361.53 


206 


50 


60 


119 



TLUSTY (Hubeny 1988) covering an effective tempera- 
ture range 10 000 < T^g < 35 000 and gravities from 
log 5 = 4.5 to the Eddington stability limit at each 
Teff. Increments were 2500 K and 0.25 dex in Toff and 



logy respectively. Models were calculated for two he- 
lium fractions, i.e. y = 0.09, and y = 0.20, where y = 
N[He]/N[H + He]. The models contain only H and He 
and assume a plane-parallel atmosphere which is in hy- 
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Table 3. Chemical abundances calculated in the model atmosphere analysis for Sher 25 and the comparison stars. The 
absolute abundance derived for each line in each star is given, and also the differential abundance of Sher 25 minus 
each comparison star. 



Differential Abundance of 
Sher 25 witli respect to 

Element/Ion Wavelength A HD2905 HD 13854 HD 14956 Sher 25 HD2905 HD 13854 HD 14956 



C II 


4267.02 


7.05 


7.27 


7.11 


7.07 


+0.02 


—0.20 


—0.04 




& 4267.27 
















C II 


6578.10 




7.12 


6.66 


6.95 




—0.17 


+0.29 


N II 


3995.00 


7.00 


7.93 


8.61 


8.19 


+1.19 


+0.26 


—0.42 


N II 


4035.08 


7.71 


8.32 


8.51 


8.40 


+0.69 


+0.08 


—0.11 


N II 


4236.98 


7.31 


7.86 


7.78 


7.84 


+0.53 


—0.02 


+0.06 


N II 


A C\ A 'i 1 O 

4241.18 


7.64 


8.07 


7.93 


8.28 


+0.64 


+0.21 


+0.35 




& 4241.78 
















N II 


A A A 

4447.03 


7.20 


7.53 


7.98 


7.45 


+0.25 


—0.08 


—0.53 


N II 


4607.16 


7.04 


8.19 


9.15 


8.91 


+1.87 


+0.72 


—0.24 


N II 


4613.87 


7.04 


8.02 


8.90 


8.45 


+1.41 


+0.43 


—0.45 


N II 


4630.54 


7.10 


8.31 


9.42 


8.53 


+1.43 


+0.22 


—0.89 


O II 


4069.62 


8.95 


8.91 


8.65 


8.39 


-0.56 


-0.52 


-0.26 




& 4069.89 
















O II 


4072.16 


9.03 


9.10 


9.03 


8.20 


—0.83 


—0.90 


—0.83 


O II 


4075.86 


8.99 


8.98 


8.88 


8.37 


-0.62 


-0.61 


-0.51 


O II 


4132.80 


9.04 


9.03 


8.93 


8.78 


-0.24 


-0.25 


-0.15 


O II 


4153.30 




9.29 


8.95 


8.95 




—0.33 


+0.00 


O II 


4317.14 


9.13 


9.06 


8.76 


8.87 


—0.26 


—0.19 


+0.00 


O II 


4319.63 


9.06 


9.12 


8.87 


8.91 


-0.15 


-0.21 


+0.04 




& 4319.93 
















O II 


4366.89 


9.22 


9.08 


8.90 


8.97 


—0.25 


—0.11 


+0.07 


O II 


4590.97 


9.28 


9.04 


9.03 


8.88 


—0.40 


—0.16 


—0.15 


O II 


4595.96 


9.19 


8.91 


9.24 


8.80 


-0.39 


-0.11 


-0.44 




& 4596.18 
















O II 


4638.86 


9.44 


9.45 


9.19 


9.37 


-0.07 


-0.08 


+0.18 


O II 


4661.63 


9.35 


9.21 


9.17 


8.75 


-0.60 


-0.46 


-0.42 


O II 


4673.74 


9.13 


8.83 


8.86 


8.80 


-0.33 


-0.03 


-0.06 


O II 


4676.24 


9.28 


9.23 


8.85 


8.95 


-0.33 


-0.28 


+0.10 


Mg II 


4481.13 


7.44 


7.58 


7.54 


7.46 


-0.02 


-0.12 


-0.08 




& 4481.33 
















Si III 


4552.62 


7.90 


8.00 


8.43 


8.17 


+0.27 


+0.17 


-0.26 


Si III 


4568.28 


7.93 


7.97 


8.50 


8.00 


+0.07 


+0.03 


-0.50 


Si III 


4575.22 


7.88 


7.87 


8.40 


7.92 


+0.04 


+0.05 


-0.48 


Si III 


4813.30 


7.35 


7.45 


7.49 


7.28 


-0.07 


-0.17 


-0.21 


Si III 


4819.72 


7.46 


7.54 


7.81 


7.54 


+0.08 


+0.00 


-0.27 


Si III 


4829.96 




7.40 


7.53 


7.37 




-0.03 


-0.16 


Si IV 


4116.10 


7.22 


7.41 


7.56 


7.43 


+0.21 


+0.02 


-0.13 



drostatic equilibrium. The consequences of these assump- 
tions have been discussed in Paper I. Line formation cal- 
culations for all the ionic species with absorption lines 
observed in the spectra of Sher 25 were calculated using 



the c odes detail (Giddings 1981), and SURFACE (Butler 
1984 ). Microturbulence was included in the calculations 



as an extra free parameter, and the statistical equilibrium 
calculations were performed for the species Hi, Hei, Cii, 
Nil, On, Sim, Si IV and Mgii. The atomic data employed 
were identical to those used in Paper I. 

Although the photospheric analysis is based on NLTE 
atmospheres, they are plane-parallel and do not account 
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Table 4. The mean of the individual hnc abundances from Table|^ for each element in each star. Also listed is 
the mean of the differential results. The absolute values are quoted on the logarithmic scale 12 + log[Nx]/[N-ii] (see 
Sect. 3.1.2). The mean abundances were determined by first determining the fractional abundance [7V2;]/[AfH] for each 
line, calculating the linear mean and quoting this on the original logarithmic scale. The errors quoted are the standard 
errors in the mean (a/y/n); where n is the number of features used, and a is the standard deviation of the sample again 
calculated on a linear scale. These represent the statistical scatter of the measurements, however a further systematic 
error will be present in each case. 



Differential Abundance of 
Shcr 25 witii respect to 



Ion 


HD2905 


HD13854 


HD 14956 


Sher 25 


HD2905 


HD 13854 


HD 14956 


Cll 


7.05 


7.20±0.07 


6.94±0.17 


7.01±0.06 


+0.02 


-0.18+0.01 


+0.16 +0.13 


Nil 


7.34±0.10 


8.09±0.07 


8.70±0.13 


8.42±0.12 


+1.29±0.16 


+0.30+0.10 


-0.13 +0.12 


On 


9.18±0.14 


9.12±0.05 


8.98±0.04 


8.87±0.07 


-0.34±0.16 


-0.25+0.05 


-0.10 +0.06 


Mg II 


7.44 


7.58 


7.54 


7.46 


-0.02 


-0.12 


-0.08 


Sim 


7.41±0.05 


7.47±0.04 


7.63±0.10 


7.41±0.07 


-0.00±0.05 


-0.06+0.05 


-0.21 +0.03 


Si IV 


7.22 


7.41 


7.56 


7.43 


+0.21 


+0.02 


-0.13 



for metal line blanketing or wind contamination. In this 
paper we concentrate on differential abundance determi- 
nations and comparing stars of similar atmospheric pa- 
rameters, and acknowledge that the absence of these phys- 
ical effects in our models may compromise the final ab- 
solute abundances. In the future it is our desire to pur- 
sue a fully unified model atmosphere and wind analysis 
of sher25 with even higher quality spectral data from, 
for example, the UVES spectrometer on the VLT. This 
would allow an excellent comparison between the very 
weak metal lines formed deep in the photosphere, and the 
stronger lines we see in the current ucles spectra. 



3.1.1. Determination of atmospheric parameters 

The relative strengths of two ionization stages of silicon 
can be used to determine the effective temperature for 
stars of early-B spectral types. Absorption lines of Sim 
and Si iv are well observed and reliably measured in the 
UCLES spectra of Sher 25. The line strengths are listed 
m Table 0, with six Imes of Si ill and one of Si iv initially 
used to determine Tegf. There is an inherent problem in de- 
tail when calculating the Si ill line strengths of the 4552A 
triplet, which has been discussed in Paper I. In the lowest 
gravity models having effective temperatures of 22 500K 
or higher, detail calculations showed an overpopulation 
of the upper levels of Sim in atmospheric regions where 
the lines are formed. This is particularly problematic for 
the 4552A triplet but is less pronounced for the 4813A 
triplet which is formed deeper in the atmosphere. The 
effect leads to a filling-in of the 4552A multiplet line pro- 
files, and hence one needs to increase the Si abundance 
substantially to offset this effect. One can readily see this 
in Table ^ were the mean abundances derived from the 
4552A triplet are between 0.5— O.Sdex higher than those 
derived from 4813A. We follow the conclusion in Paper I 
that the 4552A lines are not reliably modeled with our 



methods and use only the 4813A triplet in determining 
the effective temperature of all four stars. 

The shape of the Balmer line profiles (He, H(5 and H7, 
and H/3) are useful indicators of the stellar surface gravity. 
In particular the width of the line wings are excellent indi- 
cators of log (7. Assuming an initial Toff, the surface gravity 
is calculated by fitting the observed profiles with theoreti- 
cal lines calculated at different values of log paying par- 
ticular attention to the line wings. The mean of the results 
from H(5 and H7 was used. He suffers from interstellar con- 
tamination (Can) and a stellar blend with Hei. Both H/3 
and Ha can show significant emission contributions from 
the stellar wind and the latter is discussed in Sect. [3.2[ 
Some iteration in the simultaneous determination of Teff 
and log g was necessary, given the interdependence of each 
parameter but convergence was normally achieved within 
3-4 cycles. 

The microturbulent velocity was estimated by requir- 
ing that each line of the Sim triplet at 4552A gave the 
same abundance. The determination of Wturb in this way 
has been discussed in Paper I and it has the advantage 
that three lines from the same multiplet with quite differ- 
ent equivalent widths are used. The lines of 11 can also be 
used and these often give higher values compared to those 
derived when using Sim. Due to the homogeneity of the 
atomic data and the large range in line strengths covered, 
we use the value derived from the Si lines, and note that as 
we re-derive all parameters for the standards (and subse- 
quent abundances) in a homogeneous way any systematic 
error in this procedure will be mirrored in Sher 25 and the 
three comparison stars. 

The three parameters of Toff, logg and wturb were var- 
ied so that a simultaneous fit to all three was achieved for 
Sher 25, and the final atmospheric parameters are listed in 
Table ^ The v sin i of the star was determined by convolv- 
ing the theoretical line profiles of prominent metal lines 
with a rotational broadening function until good agree- 
ment with the observed profile was achieved (as discussed 
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in Lennon et al. 1991). We then chose three stars which 
had very similar atmospheric parameters from the list of 
bright nearby B-type supergiants analysed in Paper I. Two 
stars of almost identical Toff are HD14956 and HD13854, 
and the observed profiles of the hydrogen and helium lines 
were compared directly to make sure the luminosity class 
and the helium line strengths (and hence the atmospheric 
helium composition) were indeed very similar. This can 
be seen in Figs. |3.1.1| and p.l.l| . A third star HD2905 was 
also chosen as an extra comparison due to its relatively 
weak nitrogen line spectra, hence giving us comparison 
stars which are both qualitatively nitrogen rich and poor. 

Using the higher quality spectra now available to us (as 
described in Sect.^) we have re-derived the atmospheric 
parameters of the comparison stars using identical meth- 
ods to that used for Sher25. Our results for HD 2905 and 
HD 14956 are identical to those in Paper I, however we de- 
rive a Tes which is 1500 K cooler than the previous anal- 
ysis for HD 13854. The reason for this is likely to be the 
slightly different equivalent widths of the Si IV 4116A line 
measured in each study. We measure this line as having a 
strength of 111mA and it is quoted as 160mA in the spec- 
tra of Lennon et al. (1993). In order to keep this analysis 
as consistent as possible we use our refined temperature, 
which is similar to the previous result in Paper I within 
the quoted errors. The line profiles and strengths of the 
He I lines were all similar in Sher 25 and the 3 compar- 
isons (as seen in Fig. |3. 1.1] ). This suggests a very similar 
He abundance across the four stars within the expected 
uncertainties however there is difficulty in modeling all 
the singlet and triplet lines simultaneously and satisfacto- 
rily (see Paper I). The lines which appear to be modeled 
satisfactorily are the singlets He i 4437A and 4387A which 
form deepest in the atmosphere. However the 4437A line 
falls on top of the diffuse interstellar band at 4430A which 
is particularly strong in the Sher 25 sight line, and does not 
allow accurate extraction of the stellar line. The 4387A 
line in Sher 25 gives a helium fraction y = 0.1 ± 0.02, 
hence there is no evidence of a He overabundance. The 
447lA line, being a triplet, is difficult to model in quan- 
titative terms but the similarity in the empirical spectra 
of all four supergiants, and a measured equivalent width 
of 820A (compare with Fig. 3 in Paper I) suggests that 
there is no evidence for using a non-solar He value in the 
Sher 25 model atmospheres and that all four stars have 
comparable He abundances. 



3.1.2. Photospheric abundances 

The atmospheric parameters were used to define final 
model photospheres for all four stars. All unblended metal- 
line absorption features identified in Sher 25 (and the same 
features in the other stars) were modeled using the line 
form ation codes detail and SURFACE as described in 
Sect. |3.1.ll . A curve of growth method was used to de- 
termine an ionic abundance for each line, and the results 
are quoted on the scale 12 -I- log[Nx]/[N}i] where and 



» Hl)2905 
* Hl)13854 
□ HD14956 



0.5 - 



N Mg 

Element 



Fig. 3. Differential abundance of Sher 25 with respect to 
HD2905, HD14956 and HD13854. The values and errors 
are taken from Table ^, with the minimum error assumed 
to be 0.1 dex. 



A^H are the abundance of ion x and hydrogen by number 
(see Table H). The mean of all the individual lines were 
calculated, by first determining Nr^/Nii, calculating the 
mean and then converting back to the logarithmic scale. 
Standard deviations (cr) and standard errors in the mean 
{cr/^/n; where n is the number of features used) were cal- 
culated and are quoted in Tabled. The differential abun- 
dance of Sher 25 with respect to each comparison star was 
done on a line-by-line basis and the mean values of these 
results were calculated in a similar way. These results 
should be more reliable than the absolute abundances as 
the differential approach mitigates the uncertainties in the 
atomic data and physical methods. We discuss the results 
for each element briefly below, and the implications of 
these will be discussed in Sect.|4[ 

Carbon: The carbon abundance is based on two lines 
only, being the only reliably measured features in the 
Sher 25 spectrum. There is a Cii doublet at 6578. lOA and 
6582. 85A, however the latter is contaminated by the neb- 
ular [Nil] 6583. 6A line. There is quite strong nebulosity 
around Sher 25 which was clearly visible on the raw 2D 
CCD images, and perfect subtraction of the features was 
not possible during data reduction. Hence the 6582. 85 A 
line of Cii was not usable, although it is reliably mea- 
sured in the bright comparison stars. The 4267A dou- 
blet is well known to give unreliable absolute abundances 
with the particular line formation codes we have used (e.g. 
Vrancken et al. 2000, Lennon 1983, Paper I), hence the val- 
ues from this feature should be treated with considerable 
caution. Even in B-type dwarfs this line has historically 
proved dif ficult t o model reliably (in either LTE or NLTE) 
and Sig ut (|l996 ) suggests that the model atom of Eber & 
Butler ( 198^ ) which we have employed needs extension to 
improve results from this line. The differential abundances 
should be more reliable, and indicate that there appears 
little difference, within the errors, in the carbon abun- 



8 



Smartt et al.: The blue supergiant Sher25 






4080 4090 4100 4110 4120 
Wavelength 



40B0 4090 4100 4110 4120 
Wavelength 



4080 4090 4100 4110 4120 
Wavelength 






432B 4334 4340 4346 4352 4358 
Wavelength 



1328 4334 4340 4346 4352 4358 
Wavelength 



4328 4334 4340 4346 4352 4358 
Wavelength 



Fig. 1. The Balmer line profiles of US and H7 for Slier 25 and the three comparison stars. In each case the Sher25 
spectrum is the solid line, and the dashed profile is that from the three stars as labeled in each panel. 



dance in the four stars. The absolute abundance derived 
in all four stars is very low compared to values derived 
i n sol ar neighbourhood B-typ e star s e.g. 8.2 from Kilian 
( |l992| ) and Gies & Lambert ( 1992 ) which suggests that 
our absolute values cannot be relied upon to be used as a 
quantitative test of evolutionary models (see Sect.^. The 



only other line in common with all stars (6578A) gives sim- 
ilar results to 4267A, but Paper I found that this was only 
satisfactorily modeled at temperatures below 20000K, and 
that the hotter models have this line in emission. Hence 
with our current codes, we do not have a good measure 
of the absolute carbon abundance in these stars, although 
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Fig. 2. The line profiles of Hei 4471 A and He I 4387A for Sher25 and the three comparison stars. In each case the 
Sher 25 spectrum is the solid line, and the dashed profile is that from the three stars as labeled in each panel. 



we can say that the abundances appear fairly similar in 
each atmosphere. 

Nitrogen: Eight lines of N ii were identified and mea- 
sured in Sher 25 and each were found in the comparison 

stars, hence giving a good measurement of the nitrogen 
abundance in statistical terms. As expected, the nitrogen 



abundance in HD2905 is much lower than in either of the 
other two comparison stars and in Sher 25. This star has 
very weak N II absorption features and has been classified 
as a BCO.TIa (Lcnnon ct al. 1992). This star was chosen as 
a comparison to show the difference between the nearby 
B-type supergiants that is observed. The nitrogen abun- 
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Table 5. The errors listed in this table illustrate how 
the abundances in Sher 25 would vary within the error 
bars of the atmospheric parameters. The values were cal- 
culated around the model atmosphere with log 5 — 2.6, 
Te// = 22 300, ^ = 15kms"^ The final column lists the 
total error estimated by adding the individual errors in 
quadrature. 



supergiants. The method used to determine the mass-loss 
rate for Sher 25 is described in detail in Paper II. Exactly 
the same code is used in this case (Santolaya-Rey et al. 



Species 


logff ±0.2 


Te// ±1000 


C±5 


AlogE/H 


C II 


±0.03 


±0.02 


±0.06 


±0.07 


N II 


±0.18 


±0.10 


±0.27 


±0.34 


II 


±0.13 


±0.15 


±0.18 


±0.18 


Mg II 


±0.27 


±0.15 


±0.11 


±0.33 


Si III 


±0.18 


±0.10 


±0.10 


±0.23 



dance in Sher 25 is much greater than in IID2905, and 
is intermediate between the other two stars (illustrated 
in Fig 3). Paper I indicates that at the Teg of our stars, 
the N II line strengths in the bulk of the stars are reason- 
ably well modeled with a "normal" nitrogen abundance 
of 7.7 dex. Problems with overpopulation of upper atomic 
levels and emission features appearing in line profiles were 
not in evidence for N 11 in this temperature regime, unlike 
Cii and Sim. Hence we have more confidence in using 
the absolute values of Nii as being useful for quantita- 
tive comparisons with evolutionary calculations discussed 
in Sect.|. 

Oxygen, magnesium and silicon: With fourteen lines 
measured in Sher 25 oxygen is the best sampled element, 
as is normal in early B-type stars. Paper I found that the 
O II absolute abundances in this temperature range did 
not suffer any obvious problems with high photoioniza- 
tion rates. Given the good statistical sample of On fea- 
tures, and the encouraging analysis in Paper I we assume 
that the oxygen abundance derived in our stars is reliable 
enough to use in comparison with evolutionary calcula- 
tions. There is some evidence that the oxygen abundance 
in Sher 25 is somewhat lower than in the three standard 
stars, in particular in comparison with HD2905 and is dis- 
cussed in Section^ below where we consider the implica- 
tions of CNO processing in these objects. Both magne- 
sium and silicon in Sher 25 are very similar to the abun- 
dances derived in the standard stars. As these two ele- 
ments should not be affected by any mixing or dredge-up 
processes that may affect the CNO surface composition, 
the fact that they are similar in all three stars points to- 
wards them all having fairly similar initial chemical com- 
positions. 



3.2. Stellar wind analysis 

As the origin of the gas surrounding Sher 25 is a key fac- 
tor in the evolutionary history of this object, an investi- 
gation of the stellar wind properties is highly desirable. 
Further we can compare results directly with the exten- 
sive study of Paper II (Kudritzki et al. 1999) which in- 
vestigated the wind properties of 14 Galactic early B-type 



1997) and here we briefly outline the procedure used. 




6550 655: 



6565 
Wfivelengt.h 



Fig. 4. The Ha profiles of Sher 25 and the three compar- 
ison stars. The secondary emission peak bluewards of the 
rest wavelength of Ha in Sher 25 is due to imperfect neb- 
ular subtraction in the 2D image. Weak residuals of the 
[Nil] 6583. 6 A & 6548.1 A nebular lines can also be seen. 
The Ha emission of Sher 25 is similar in strength to the 
comparison stars and the wind analysis implies it has quite 
normal mass-loss 

The stellar parameters determined above were used to 
calculate the emergent flux from a non-LTE model atmo- 
sphere, and by assuming an absolute visual magnitude of 
the star, a radius can be determined. However the cal- 
culation of absolute V-band magnitude depends on our 
knowledge of distance and reddening. The former has been 
calculated by several auth ors, th e most recent of which is 
6.3±0.6kpc (Pandeyet al. |2000|) . Sher 25 has B-F = 1.4, 
which indicates a extinction E(B-V|=1.6 from the intrin- 
sic colours of Deutschman et al. ( |l976 ). Pandey et al. also 
determine R (the ratio of total to selective reddening) to- 
wards NGC3603 and suggest a value of 4.3, which is much 
higher t han th e normal Galactic line of sight value R=3.1 
(Seaton 1979 ). However this may be applicable more to- 
ward the central regions of the cluster, whereas Sher 25 
sits outside the main body of young 0-stars. We experi- 
mented with the value of 4.3 and determined what effect 
it would have on the derived intrinsic luminosity of the 
star and its radius and whether this was consistent with 
the observed Ha profile. Assuming R=4.3 would imply 
a luminosity of log L/Lq — 6.2, and a radius of 8OR0. 
These values suggest the star would be more luminous 
and have a larger radius than HD190603 (Paper I, Lennon 
et al. 1992), but the surface gravity of Sher 25 we derive 
is incompatible with a very low gravity, extremely high 
luminosity star. HD190603 is probably the most luminous 
star in the Paper I sample, and shows Ha very strongly 
in emission. Sher 25 has weaker Ha emission and broader 
absorption wings in the other Balmer lines, which strongly 
suggests it is a lower luminosity object. Also a luminosity 
of 6.2 would suggest that Sher 25 is much more luminous. 
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and hence has a much lower surface gravity than either 
of HD2905 or HD13854, both of which have undergone a 
wind analysis in Paper II. However the Ha emission line 
in Sher 25 is not significantly stronger than in either of 
the latter two, and the Balmer line profiles do not suggest 
a much lower surface gravity. Adopting a value of R=4.3 
would lead to serious inconsistencies between the luminos- 
ity of Sher 25 and its spectral properties, hence it is very 
unlikely that such a high value of R is applicable. There 
does appear to be evidence for a variation in R toward this 
sight line, and we adopt a value of R=3.7 as being indica- 
tive of the true value with a probable error of ±0.5. With 
this value we obtain a radius of 59 ± 15i?0. The observed 
Ha profiles of the stars are shown in Fig.Q. The emission 
peak on the blue wing of Sher 25 is indicative of imperfect 
nebular subtraction, however the rest of the profile will not 
be affected given its width and the narrow nature of the 
nebular line. Crowther & Dessart ( 1998 ) have also derived 
a distance to NGC3603 based on the absolute magnitudes 
of 0-stars. They find a distance of 10.1 kpc, using the stan- 
dard value R=3.1. We find that Using these numbers puts 
luminosity and radius of Sher 25 at very similar values to 
those derived using d=6.3kpc and R=3.7 i.e. 59 ± 15i?0 
and logL/Lo = 5.9 ± 0.2. 

No UV spectra exist to determine a terminal veloc- 
ity for the wind of Sher 25 we adopt a value of 1000km 
s"^ The values of v^o for HD2905 and HD13854 are de- 
rived in Paper II and values of 1100 and 1000km s~^were 
found, indicating our adopted value for Sher 25 is reason- 
able. The mass-loss rate M can be obtained from fitting 
the Ha profile, using the non-LTE unified model atmo- 
sphere technique incorporating the stellar wind and spher- 
ical extension, as described in Paper II. A grid of models 
was calculated with different mass-loss rates and the best 
fit was found for 1.75 x lO~^M0yr~^, and /? = 1.5 (see 
Fig.||). The mass-loss rate is close to (and roughly inter- 
mediate between) the values for HD 13854 and HD2950 in 
Paper II, and appears fairly typical for early-B supergiant 
stars with solar-type chemical compositions. The fit to 
the emission profile is good, but the absorption dip blue- 
ward of the peak is not reproduced. This dip is also seen 
in other early B-types of similar atmospheric parameters 
and mass-loss rates, and the failure of the wind analysis 
calculations to match the observed shape is discussed in 
Paper II. It may be due to the neglect of metal-line blan- 
keting in the models, or wind variations and deviations 
from spherical symmetry. We further derive a luminosity 
of log L/Lq — 5.9 ± 0.2, and a modified wind momentum 
of log(Muoo(i?/i?o)"'^) = 28.9 ± 0.3, which would place 
the star close to the wind momentum-luminosity relation- 
ship for solar metallicity early B-type stars (see Paper II). 
Hence the most important conclusion from the fitting of 
the Ha and derivation of its wind and luminosity proper- 
ties is that the wind of Sher 25 looks fairly typical for an 
early B-type supergiant. 
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Fig. 5. The Ha profile of Sher 25 and the fit from the 
unified wind models. The peak is well reproduced, whereas 
the bluewards absorption is not seen in our model spectra 



due to the reasons described in Sect. 3.2 



4. Discussion of the evolutionary scenarios for 
Sher 25 

4.1. A post red-supergiant star in a "blue-loop" 

The existence of the ring shaped nebula and bi-polar out- 
flow around Sher 25 has been interpreted by Brandner et 



al. (1997a, 1997b) as arising during a previous red super- 
giant evolutionary phase. In this scenario the star began 
its life as an 0-type main-sequence object, burning H to 
He through the CN and the ON-cycles. After core hydro- 
gen exhaustion, this evolutionary picture involves the star 
moving quickly across the HR diagram to the red super- 
giant region where it undergoes core He burning. During 
the red supergiant phase the stellar envelope becomes fully 
convective, which would mix a certain amount of the core 
material into the envelope. The relative abundances of 
CNO in the red supergiant atmosphere would then change 
drastically. The predicted final ratios of these elements will 
clearly depend on the mass of CN and ON-cycled core 
gas which is mixed into the envelope (assumed to have 
the stars natal composition), and the mass of the enve- 
lope which is left at the beginning of the red supergiant 
phase. The envelope mass will be affected by the amount of 
mass the star has lost during its hot main-sequence phase 
through radiatively driven winds. This global picture is 
simply and clearly described in Lamers et al. (2001). A 
single major mass-loss event is then invoked to explain 
the nebula and this dense, cool gas should have a chemi- 
cal composition similar to the polluted red-supergiant at- 
mosphere. Hence the abundances in the ionized nebula 
around Sher 25 could be used to test the validity of this 
argument, however there has been no accurate determina- 
tion of abundances in the nebula to date. The spectra pre- 
sented by Brandner et al. (1997a and 1997b) does not have 
enough wavelength coverage to estimate electron densi- 
ties and temperatures consistently, hence no quantitative 
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Table 6. Comparison of CNO abundances from our analysis of the four B-type supergiants. For comparison the values 
for solar neighbourhood B-type main-sequence stars are taken from Gies & Lambert (1992), and the solar values are 
from Grevesse & Sauval 199^ The C, N and O element abundances are by number fraction. The numbers in italics 



should be treated with extreme caution, given the uncertainty in the absolute C abundance as discussed in Sect. 3. 1.2 



Element Solar B-stars HD2905 HD 13854 HD 14956 Sher25 



c 

N 


C+N+O 
H + He 


0.33 
0.08 
0.68 
0.99 


Element abundance by number (xlO 
0.16 0.01 0.02 0.01 
0.07 0.02 0.12 0.50 
0.48 1.51 1.32 0.95 
0.65 1.41 1.33 1.33 


0.01 
0.26 
0.74 
0.92 


N/C 
N/0 


0.25 
0.12 


0.41 
0.14 


Elemental fractions 
2.0 7.7 
0.02 0.09 


55.7 
0.53 


26.3 
0.36 



nebular abundances can be derived. These spectra do in- 
dicate that the [Nil] /Ha line strength ratio is somewhat 
higher than the background H ii region, which Brandner 
et al. interpret as a sign of real nitrogen enhancement in 
the circumstellar gas, however a full and accurate nebular 
analysis is still required. After a short period as a red su- 
pergiant, and after the mass- loss has occurred it has been 
speculated that Sher25 (and other blue supergiants) re- 
turn to the blue region of the HR diagram, undergoing 
what have been called "blue-loops" . 
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Fig. 6. Stellar evolutionary tracks for stars with zero-age 
main sequence masses of 20, 25, 40, 60 and 85M0from 



Meynet et al. |1994 The symbols are Sher25 (filled star), 
Sk-69°202 progenitor of SN1987A (filled circle), HD2905 
(open diamond), HD13854 (open square), HD14956 (open 
triangle). Luminosities are taken from Paper II and 
Walborn et al. (1989); apart from IID13854 which was 
calculated assuming membership of Per OBI, distance as 
in Paper II, and optical magnitudes from the Simbad 
database. 



We can use the photospheric abundances we have de- 
rived for Sher25 to test this model of a mass-loss event 



occurring when the star was a red supergiant and sub- 
sequent bluewards evolution. Lamers et al. (2001) have 
recently presented comparisons between chemical abun- 
dances in LBV nebulae and what one would expect in 
different evolutionary scenarios. They come to the conclu- 
sion that LBV nebulae are ejected during the blue super- 
giant phase, and that such stars have not gone through 
a previous red supergiant episode. We follow their line 
of argument to determine if the CNO abundances we de- 
rive in Sher 25 are consistent with a previous red super- 
giant phase. The Si and Mg abundances in Sher 25 are 
very similar to those in the three solar neighbourhood 
comparison stars, and these are not expected to be af- 
fected by any mixing of H-burning products into the pho- 
tosphere. The Galactocentric distance of NGC3603 is ap- 
proximately 8.5 kpc, roughly the same as the Sun, and 
galactic abundance gradients derived in RoUeston et al. 
( ^OOOD and Smartt & RoUeston ( |1997| ) indicate that there 
is fairly low scatter in stellar metallicity at this distance 
from the Galactic centre. Hence we assume that Sher 25 
had an initial chemical composition similar to B-type stars 
and nebulae in the solar neighbourhood. Table|| lists the 
number fractions of CNO derived in Sher 25, the three 
comparison supergiants, the Sun, and solar neighbourhood 
B-type stars. From Table^ and Table ^ it is apparent that 
Sher 25 has a nitrogen abundance approximately a factor 
of two lower than HD 14956, and a factor of two higher 
than HD13854. These two stars have been classified in 
Paper I as "highly processed" in that they are amongst 
the group of Galactic supergiants that show the highest 
N/C line strength ratios. The line strengths of Sher 25 
are comparable to stars in this category, and the abun- 
dances we derive place it intermediate between HD 14956 
and HD 13854 in terms of N enrichment. However it is 
certainly not extreme in its nitrogen abundance. It has a 
much higher abundance than HD2905, however we chose 
this star as a comparison specifically because it was one of 
the objects in Paper I with the lowest apparent N enrich- 
ment. This illustrates the wide variation in nitrogen abun- 
dances that is seen amongst early B-type supergiants, and 
the abundance in Sher 25 is similar to those Galactic stars 
which appear most N enhanced. 
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Lamers et al. (2001) have presented the envelope abun- 
dances that are predicted by the models of Meynet et al. 
(1994) after convective mixing in the red supergiant phase. 
This depends critically on two factors. Firstly it depends 
on the mass- loss during earlier phases, which is dominated 
by radiatively driven winds while the star is hot. Clearly 
the higher the mass-loss, the more of the envelope is lost 
and hence there is less gas of initial composition to mix 
with the CNO-processed core material. Secondly it will de- 
pend on how far into the core the convection reaches and 
how much core gas actually gets dredged up into the con- 
vective atmosphere. Lamers et al. (2001) have plotted the 
N/0 and N/C ratios for various values of mass-loss rates 
and initial stellar mass. We estimate the zero- age main- 
sequence mass of Sher 25 to have been 6OM0 from com- 
parison with the Meynet et al. tracks in Fig. |[ Although 
we have directly derived a mass-loss rate for Sher 25 in 
its present status as a luminous B-type supergiant, it may 
have had a higher M when it was a main sequence star. 
A 60Mq star on the main-sequence would have a mass- 
loss rate of typically ~ 2 — 5 x lO^^M© yr~^ (Puis et 
al. 1996), hence during its H-burning lifetime (4Myrs; 
from the Geneva models) it would have lost 8-2OM0 of 
its envelope through its stellar wind. The surface N/0 
and N/C ratios show a steep dependence on M in this 
range, and the Lamers et al. calculations indicate ranges 
of 3 ^ N/0 ^ 20 and 6 ^ N/C 30. 

In order to compare this meaningfully with the abun- 
dance ratios we derive, we must have confidence in the 
absolute abundances derived in our methods. To date the 
evaluation of this has been difhcult, as pointed out in 



Smartt et al. (2001a, 1997) and Paper I. These studies have 
concentrated on differential analyses and comparisons of 
the difference between stellar subgroups, rather than as- 
suming validity of the absolute abundances. In particular 
the absolute value of the carbon abundance in our stars is 



of dubious validity, as we discussed in detail in Sect. 3.1.2 . 
Hence we will not quantitatively compare the N/C ratios 
in our stars to those in the evolutionary models. However 
as discussed in Sect. 3.1.2 the N and O absolute results 



should be reliable enough to carry out a meaningful com- 
parison, providing we consider the probable errors in the 
methods. In Table ^ the N/0 fraction is listed for each 
star. Given the uncertainties in the individual abundances, 
these ratios are uncertain within a factor of approximately 
3. For Sher 25, the observed N/0 ratio (0.36) is much lower 
than that expected from the Lamers et al. (2001) calcu- 
lation for a star which has undergone convective dredge 
in the red supergiant phase (3—20), even allowing for an 
uncertainty of 3 in our derived results. For the most ni- 
trogen rich star in our sample HD 14856 (of initial mass 
approximately 4OM0) N/0 is predicted to be greater than 
2, compared with our estimate of 0.53. This also implies 
the observed N/0 ratio is too low to be compatible with 
convective dredge-up having occurred. The global analy- 
sis of the large sample of B-type supergiants in Paper I 
implied that their CNO line ratios and estimated abun- 
dances were not consistent with them (or even a subset 



of them) being post red supergiant stars. Our analysis of 
Sher 25 indicates that it has similar CNO abundances to 
the "highly processed" stars, but that it is not extremely 
nitrogen rich. This suggests that, despite the appearance 
of the nebula, the star has not undergone dredge-up and 
thus has not previously been a red supergiant. 

4.2. Evolution direct from tine main sequence 

Although we see no strong evidence for a previous red- 
supergiant stage of Sher 25, the abundances do suggest 
that some form of CN and NO cycle processed gas has 
been mixed into the stellar atmosphere. The sample in 
Paper I indicates that there is a range in the CNO abun- 
dances in stellar atmospheres. This is consistent with rota- 
tionally induced mixing while the stars are main-sequence 
0-types and/or as they evolve off the main-sequence to the 
cooler surface temperatures that we presently see. In this 
case one would envisage a range of C/N and N/0 ratios 
in the atmospheres due to the distribution of rotational 
velocities. Lamers et al. (2001) have also estimated the 
C/N and N/0 surface abundances for massive stars that 
undergo rotationally induced mixing. They calculate the 
surface abundance ratios as a function of the mixing time 
Tmix and find that their results agree well with the calcu- 
lations from Heger & Langer ( |2000| ) and Maedcr ( |1998D if 
Tmix — 1 — 3rMs (where tms is the main-sequence lifetime). 
Further they estimate that Tmix = 2-4 tms is needed to ac- 
count for their observed LBV nebular abundance ratios. 

For a 6OM0 star the N/0 ratio in the atmosphere of a 
blue supergiant which has undergone rotationally induced 
mixing is in the range of 1 <N/0< 3 if Tmix is 2 — 4 
times the main sequence lifetime. The value for Sher25 
(N/O=0.36) is lower than these predicted values for 2 < 
Tmix < 4, but is a factor of 3 higher than the expected 
initial value of N/O==0.12. Our measured value suggests 
that Sher25 has undergone less mixing than these models 
predict, although almost certainly some moderate amount 
of mixing has occurred to enhance the N abundance. This 
suggests a larger mixing time (i.e. greater than 2 — 4) is 
more appropriate for Sher25. 

5. Conclusions and wider implications 

5.1. Implications for the evolutionary status of blue 
supergiants 

From our detailed abundance analysis of Sher 25 and com- 
parison to evolutionary calculations we can find no clear 
evidence that the star was previously in a red super- 
giant phase. The CNO abundances we derive at its surface 
are similar to those in other Galactic B-type supergiants 
which show somewhat enhanced N/C and N/0 abundance 
ratios. However the extent of this enhancement does not 
appear to be high enough to suggest the star dredged up 
CN and NO cycle equilibrium core gas during a cool su- 
pergiant phase. Even allowing for a large mass-loss rate 
during its main-sequence lifetime does not reconcile the 
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results, and the current wind properties of Slier 25 also in- 
dicate that it is a relatively "normal" B-type supergiant. 
Certainly this star sits comfortably in the largest sam- 
ple of Galactic B-type supergiants yet analysed (Paper I 
and Paper II), in terms of its mass, atmospheric parame- 
ters, abundances and wind properties. It does appear to 
be similar to the group of supergiants which show the 
largest N/C line strengths, and hence N/C abundance ra- 
tios but our quantitative comparison of the N/0 results 
are more compatible with it having undergone some form 
of rotationally induced mixing either during or just after 
its main-sequence phase. 

The sample of Galactic A-type supergiants analysed 
by Venn ( 1995 ) shows evidence for the stars having under- 
gone partial mixing of CN-cycled gas. However the abun- 
dances derived are not suggestive of the stars having been 
through a first dredge up and hence are unlikely to be 
post red supergiant objects. Venn concludes that these 5- 
20Mq stars have evolved directly from the main-sequence. 
The Galactic B-type supergiants analysed in Paper I are 
of higher mass (~20-40Mo), and the results are also not 
consistent with them having undergone a previous red su- 
pergiant evolutionary phase. The study of the abundances 
in the nebulae around LBVs (in the 40-80Mq mass range) 
by Lamers et al. (2001) again suggests that they have not 
previously been red supergiants and that the nebulae are 
ejected during the blue supergiant phase. The chemical 
enhancements seen in these three cases are possibly all 
due to some form of rotation induced mixing, which will 
have some broad distribution based on initial rotational 
velocities and mass-loss rates. 

We find that the blue supergiant Sher25 is also very 
unlikely to have gone through a red supergiant phase. 
Compared to other B-type supergiants in the galaxy, its 
atmospheric chemistry and wind properties are not pe- 
culiar in any way. Its peculiarity stems from the exis- 
tence of the striking circumstellar nebula. Lamers et al. 
(2001) have speculated that the nebulae around LBVs 
have been ejected in the blue supergiant phase, although 
what the mechanism actually responsible for this is open 
to some debate. Therefore in the Galaxy, there is no sin- 
gle blue supergiant or nebula surrounding a blue super- 
giant, which displays abundances consistent with the star 
having evolved through a previous red-supergiant phase. 
The idea that nebulae surrounding massive, evolved stars 
must have been ejected during the red supergiant phase 
does not appear to be supported by stellar analyses of 
blue supergiants shown here, in Paper I and in Lamers et 
al. (2001). 

5.2. Comparison with SN1987A 

The origin of Sher25's nebulae is thus not well under- 
stood. Brandner et al. (1997b) have estimated the dynam- 
ical age of both the equatorial ring and bi-polar lobes to 
be around 6600 yr, and the total ionized gas to be be- 
tween 0.3 — 0.7Mq. The current mass-loss rate is much 



too low to have had much impact on the mass of this 
material, as in 6600 yr the mass of gas ejected in the 
current fast wind is ~O.O1M0. Brandner et al. (1997a, 



1997b) have also speculated that the Sher25 nebula is 



similar in morphology, mass, composition and kinematics 
to that now seen around SN1987A. The SN1987A nebulae 
was almost certainly there before explosion, obviously too 
faint and compact to have been detected on pre-explosion 
plates, suggesting the progenitor star Sk— 69°202 was at a 
similar evolutionary state to Sher 25 (although somewhat 
lower mass). By implication our results then suggest that 
Sk— 69°202 may never have gone through the red super- 
giant phase, but rather evolved directly from the main- 
sequence and ejected its circumstellar material during the 
blue supergiant phase. 

The initia l ejecta of SN1987A was studied by Fransson 
et al. (|1989|) with lUE and ratios of N/C~7.8 and 
N/0~1.6 were found. The latter is significantly more than 
that in the atmosphere of Sher 25 and could be compatible 
with the progenitor having had a previous red supergiant 
phase from comparisons with the Lamers et al. calcula- 
tions. However the modeling of this emission line spec- 
trum is problematic and uncertainties in N/0 are proba- 
bly of the order ±1. Three features of the SN1987A event 
have caused problems for theoretical models looking for 
a consistent explanation of the event; the relatively high 
temperature of Sk— 69°202; the stellar ring around the SN 
(Wampler et al. 199C, and Jakobsen et al. |l99l| ); and the 
chemical anomalies discussed by Fransson et al. (1989). 
The review by Podsiadlowski ( 1992| ) suggests that no cur- 
rent theory of the evolution of single stars can simulta- 
neously cope with all three and that binary accretion (or 
merger model) is theoretically the most promising, espe- 
cially with regard to the mechanism for ring formation. We 
find no evidence for radial velocity variations of Sher 25 
(within the uncertainties of ~7km s~^) although the four 
1200s exposures were taken back to back. We plan future 
velocity monitoring of this star to determine whether or 
not it is a binary system. In the wider picture the ques- 
tion of how many blue supergiants are supernova Type II 
progenitors is very uncertain although new projects with 
data mining and the rich data archives which now exist 
may shed some light on this in the future (e.g. Smartt et 
al. |2001b ). It would be also worthwhile to determine the 
CNO abundance ratios in other early B-type supergiants 
in the LMC in order to gain insight into the stars of similar 
metallicity at a similar evolutionary state to Sk— 69°202. 
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